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REPORT  OF  PROGRESS  IN  THE  DEVELOPMENT  OF  AN 
ELECTRON  IMPACT  IONIZATION  MODEL 

I . Introduction 

The  radiative  emissions  from  high  temperature  laboratory  and  solar 
plasmas  contain  a potential  wealth  of  information  about  its  chemical 
composition,  temperature  and  density  structure,  and  radiative  cooling 
rate.  However,  before  any  reliable  quantitative  information  can  be 
deduced  from  the  plasmas'  spectral  emission  features,  some  knowledge 
must  be  available  on  the  fundamental  processes  leading  to  the  spectrum 
formation.  In  principle  this  involves  the  construction  of  a plasma 
model  that  includes  all  the  relevant  and  competing  physical  processes 
necessary  to  describe  the  dynamic  evolution  of  the  spectral  emission 
features.  Prominent  among  these  processes  is  electron  impact  ionization 
and  its  inverse  process,  three-body  recombination. 

The  theoretical  study  of  electron  impact  ionization  of  atoms  and 
ions  has  been  under  intensive  investigation  for  a period  of  years.  How- 
ever, only  in  recent  years  has  substantial  progress  been  achieved  with 
the  aid  of  large  computers.  Nevertheless,  the  majority  of  calculations 
employ  either  the  Born  approximation  or  some  variation  thereof  and  a 
variety  of  semi-classical  and  empirical  methods.  It  is  well  known  that 
outside  the  threshold  region  the  accuracy  and  validity  of  the  Born 
approximation  becomes  better  the  higher  the  incident  energy.  Unfortu- 
nately, the  magnitude  of  the  ionization  rate  is  often  dominated  by  the 
threshold  region.  Hence,  the  development  of  semi-classical  and  empirical 
methods  to  attempt  to  characterize  the  threshold  region.  Of  the  many 

methods  available  the  Exchange  Classical  Impact  Parameter  (ECIP)  method 
Note:  Manuscript  submitted  February  10,  1978. 
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has  met  with  some  success,  but  it  is  not  yet  fully  understood  why.  For 
charge  states  greater  than  about  25  Lotz's  empirical  results  are  often 
applied. 

The  purpose  of  this  report  is  to  summarize  our  progress  to  date  on 
the  development  of  an  electron  impact  ionization  formalism  and  computer 
code.  The  collision  process  is  treated  quantum  mechanically  using  the 
method  of  distorted  waves.  The  results  of  some  benchmark  calculations 
with  hydrogen  will  be  presented. 

II.  Outline  of  the  General  Theory  of  Electron  Impact  Ionization  of  Atoms 
We  shall  consider  the  ionization  of  an  atom  or  ion  with  n equivalent 
electrons  and  we  shall  assume  that  the  conditions  of  the  LS  coupling  are 
satisfied  for  all  atomic  electrons.  Atomic  units  (e  = m = *.)  will  be 
used  throughout  this  report. 

The  total  cross  section  Q(£nSL;k  ) for  ionization  of  an  atom  in  the 

o 

initial  state  £nSL  by  collisions  with  electrons  having  incident  energy 

. 2 

\ kQ  and  random  spin  orientation  is  given  by 

QanSL;k  ) = Z q UnSL;k  £n“1S'L')  (in  atomic  units  a2)  , (1) 

o S,L,  o o 

where  q is  the  cross  section  for  the  ionization  processes  in  which  the 
final  state  of  the  ion  is  specified  by  quantum  numbers  S'L'  . After  an 

-4  -4 

ionizing  collision,  the  two  free  electrons  have  momenta  k^  and  k„  and 
because  they  cannot  be  distinguished,  it  is  not  possible  to  say  which 
one  is  the  scattered  and  which  one  is  the  ejected  electron.  Without 
the  loss  of  generality  we  therefore  adopt  the  convention  k3  S k^.  In 
analogy  to  the  classical  description  of  ionization,  the  electron  having 
momentum  kQ  will  be  called  henceforth  the  "ejected"  electron. 


The  cross  section  q can  be  written  as  an  integral  over  the  energy 
of  the  ejected  electrons'', 

E/2 

qanSL;k  = k_1  f k.  k^cr(  2nSL  ; k -*£n"1S'L'  ;k.  k.)  d kh  , (2) 

O O J l d O Id  d 

O 

p 

where  E = 4 k - E.,  E.  being  the  ionization  energy,  and  g*  is  given  by 
e o 11 

the  scattering  amplitudes  f describing  the  individual  ionization  pro- 
cesses, 


(KAUSL*!c  -*£n’1S'L'  ;k1k  ) = [3iK2S+l)  (21+1)'1 


-1 


M^u 


V*L  v 


dk  dk.  dk_. 


(3) 


In  (3),  k , k.  and  k~  are  unit  vectors  in  the  direction  of  k , k, , k. 

’ o 1 2 o 1 2 

and  y,  , n, , ■ ■ specify  the  spin  orientations  of  the  free  electrons 
0 12 

before  and  after  the  coll is  ion , ] f 1 ^ is  integrated  over  all  orientations 
of  k^,  ic^ , summed  over  all  final  states  characterized  by  MS  «L  “lu2 
and  averaged  over  all  directions  of  incident  electrons  and  all  initial 


states 


*s\*< 


The  scattering  amplitudes  f are  defined  by  the  asymptotic  form  of 

the  total  wave  function  Y describing  the  scattering  of  electrons  with 

2 

incident  energy  \ kQ  . The  asymptotic  form  of  w for  the  case  of  two 

2 3 

electrons  in  the  field  of  an  ion  was  derived  by  Peterkop  ’ and  by 
Rudge  and  Seaton4-.  The  knowledge  of  the  correct  form  of  in  the 
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asymptotic  region  permits  us  to  write  an  explicit  expression  for  the 
scattering  amplitude  as  an  integral 


ko  - kjk2> 


= . eiA(^lk2)  n Y (xg xn+]>;  x^CHd.S)  + E'-E] 


xCY*(/"1s'L'Ms'HL';  ?3— -xn+i)  fv(^;  0'1 ) Xj_  cz1-i?1;  ?x)  (k) 


Xv(  Ug’  *2  ^2_k2’  r2^  ” v c U- L ; Cfp)  X]_  (Z^-k^;  r^)  y (^2> 


X x2(Z2-k2;  ?1)]]  d*l d*n+l 


H(l,2)  is  the  hamiltonian  of  electrons  labeled  1,2,  E'  the  energy  of  the 

final  state  of  the  ion,  E the  total  (positive)  energy  of  the  system  and 

x.  represents  spatial  and  spin  coordinates  of  the  i-th  electron. 

v (x„  ....x  . ; x, ) is  the  exact  solution  of  the  Schrodinger  equation 

— 2 n+L  1 

HY  = EY  (5) 


and  it  is  normalized  so  that  for  r^  -*  ® it  has  the  form 

I - - \ — Vl>  v <vV  tw2,.^(Vi-Vi>rl1 

x exp  [ilc^-i  (Z/kQ)  log  (k^-ic^  )] 


(6) 


+ scattered  waves 


for  collisions  with  atoms  or  ions  with  asymptotic  charge  Zv.7=0  for 


neutral  atoms),  w is  antisymmetric  with  respect  to  x?  ....  only. 


(u^CTjP  in  (9-),  (5)  and  (6)  is  a spin  function  representing  spin 


orientation  u,  and  cr^  is  the  spin  coordinate  of  the  k-th  electron,  w 


is  the  antisymmetric  function  of  the  atom  before  or  after  the  ionizing 
collision. 


The  functions  XjCZj-kjjr)  satisfy  equations  of  the  form 


(7) 


where  V^(r)  may  be  any  short-range  potentials.  At  large  r. 


X .(Z  -lc.  ;r)~ri+Z^fik^(k  .r+ic  .r) } ^ exp  r-i£.r-i(Z  ./k .)  log  (k,r4k.r)1 

i i : ~ r j j j J J j j j j j 


(8) 


The  charges  Z^  and  Z^  in  (k) , (7)  and  (8)  are  related  to  the  phase  factor 

-4  -4  -4 

exp  [i&(k^kg)]  in  (9-)  which  is  divergent  unless  Z^,  Z , k^  and  k„  satisfy 
the  equation 


Z1  Z2  _ Z+l  , Z+l  1 

kl  k2  “ kl  k2  |^-k2! 


(9) 


In  (4)  we  have  omitted  the  term  representing  "classical  exchange", 
i.e.  a capture  of  the  colliding  electron  and  a subsequent  ejection  of  two 
atomic  electrons.  This  process  has  been  described  in  Ref.  1. 

Eq.  (k)  is  the  basis  for  the  calculation  of  scattering  amplitudes  f. 
Cross  sections  can  be  obtained  using  relations  (1),  (2)  and  Q).  However 


the  exact  solution  of  the  Schrodinger  equation  (5)  is  unknown,  and 
therefore  y has  to  be  replaced  by  an  approximate  function.  This  pro- 
cedure leads  to  various  approximations.  Furthermore,  if  y is  not  the 
exact  solution,  f becomes  dependent  on  the  choice  of  potentials  V and 
Vg  in  (7).  In  actual  calculation  it  is  difficult  to  satisfy  eq.  (Q)  for 
effective  charges  Z^  and  Z^  which  depend  on  relative  directions  of  k^ 
and  kg.  If  the  requirement  (9)  is  not  satisfied,  the  phase  factor 
exp  riA(k^kg)]  is  essentially  undefined  and  then  there  exists  an  arbi- 
trariness in  the  relative  phase  assigned  to  the  two  terms  in  ” — - j in 
eq.  (k) , the  first  giving  the  direct  part,  the  second  the  exchange  part 
of  the  scattering  amplitude. 


Ill . Partial  Wave  Theory  and  Approximations 

For  numerical  calculations  it  is  convenient  to  expand  the  function 
y.(Zj-lCj;r)  into  partial  waves  using  the  relation'' 


1 . 03 

(kj)  * r"  2T  \ exp  [i(T^  + Fk  z (r) 

i—O  j 


(10) 


hi 

m -- 1 


* E <-G.)  (9)  , 


where  are  spherical  harmonics  defined  by  directions  of  -k^  and 

r,  respectively,  T|  is  the  Coulomb  phase  T|  = arg  T (l  + 1 - iZ,/k.) 

I l J J 

and  F,  is  the  solution  of  the  equation 

r 
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± 


Ai 


r _d^  1(1  + 1) 

- 2 ‘ 2 

dr  r 


with  the  asymptotic  form 


2Z  . 


2-i 


+ 2 V (r)  + k,]  F (r)  = 0 
r J J Kj* 


(ID 


F,  . ~ k.'s  sin  [k.r  - \ + (Z./k.)  log  (2k  r)  + Tj  + A.].  (12) 

kjX  J J J J J * 1 


A.  is  the  phase  shift  depending  on  V.(r)  and  * = 0 if  V . = 0. 

I 3 A J 

We  adopt  the  following  approximations: 

(a)  The  wave  function  V of  the  total  system  appearing  in  (U.) 
is  replaced  by 


C;W  *c 


(Zko;ri) 


(13) 


where  ^(zi^jr^)  is  defined  similarly  to  (7)  and  (3),  and  Z is  the 
asymptotic  charge  of  the  atom  or  ion  before  the  collision.  We  assume 
that  the  radial  functions  of  atomic  electrons  remain  unchanged  by  the 
ionizing  collision.  The  potential  V (r)  in  equation  (7)  for 
Xo(Z.kQ;r^)  represents  the  modification  of  the  Coulomb  potential  Z/r 
near  the  nucleus  and  in  our  approximation  it  is  given  by 


VQ(r)  = nr'1  - n Q(r) 


with  Q(r) 


- t_1  / P£  * J 


P£ 


(1M 

(15) 


where  n is  the  number  of  atomic  electrons  before  the  ionizing  collision 
and  P (r)  is  the  radial  function  of  atomic  electrons. 

As 


1 


(b)  We  will  assume  that  charges  Z^  and  Z^,  which  determine  the 
asymptotic  forms  of  the  wave  functions  '/^  and  of  the  ejected  and 
scattered  electrons, are  not  dependent  on  the  direction  of  and  k, . 

This  assumption  substantially  simplifies  the  evaluation  of  the  cross 
sections,  but  it  is  in  violation  of  requirement  (9).  In  the  majority  of 
calculations  on  ionizing  collisions  it  has  often  been  assumed  that  Z^  = 

Z and  Zg  = Z + 1.  We  intend  to  study  the  effect  of  different  choices 
of  Z^  and  Zp  including  non- integer  values  of  these  charges.  The  short- 
range  potentials  and  in  (7)  may  be  arbitrary  and  the  result  will 
depend  on  their  form.  One  of  many  possibilities  is  to  set 


Vj  = (n-£\)  [r  1 - Q (r)  ] 


(16) 


with  Q(r)  defined  by  (15). 


such  that 


Z . 
J 


Sj 


(17) 


and  0 ^ £.5  1.  This  particular  choice  of  and  would  make  our 
approximation  similar  to  the  distorted-wave  method  in t he  calculation 
of  collisional  excitations.  Born  and  Coulomb-Born  approximations  would 
be  obtained  if  we  put  VQ(r)  = V^(r)  = V^(r)  =0,  Z^  = Z,  Zp  = Z + 1. 

(c)  We  omit  exchange  contributions  to  the  scattering  amplitude, 

i.e.  we  ignore  the  second  term  in  '* ] in  eq.  (k) . Effects  of  exchange 

will  be  studied  later.  In  a non-exchange  approximation  the  value  of  the 
phase  factor  exp  [iACk^k^)]  is  irrelevant. 

With  the  assumptions  (a),  (b)  and  (c),  the  calculation  of  the  ioni- 
zation cross  sections  is  considerably  simplified  and  the  problem  is  in 
many  respects  analogous  to  collisional  excitation. 
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The  function  X (Zk  ;r.)  may  also  be  expanded  into  partial  waves 
o o i 

according  to  (1C),  where  for  convenience  we  write  Y*  (k  ) Y,  (r,  ) 

2 m o 2 m 1 

A A 0 0 0 0 

instead  of  Y,  _ (k  ) Y*  m (r  ) . 

* tn  o 2 m L 
0 0 0 0 

Using  the  notation 


-1 


VVj^r^)  3 v (y‘j;CTi)  ri  Y^j^j(ri) 


(16) 


for  a one-electron  wave  function  with  spin  the  expression  (k)  for 
the  scattering  amplitude  takes  the  form 


,l!”SLasaL>‘o-  s=  - V2=  ei  e2 


= - (2irf  exp  * i ' 1 n 2 (W) ' (k^kg) 


i +2,  +2, 


x£  1°  1 2expr10,  «,  + * +«  +,  )]  £ (£  > 

, „ „ O i .2  o L d 00 


£o^2 


m m,  nu 
o 1 2 


(10) 


Xh1»1(-h)Yt2m2(-62)|*«"siMsML; 


xo x„,i)  (2  m u k ;x  ) 

d.  n+J.  O 0 0 0 0 i 


X[H(1,2)+E’-E]  Y*(An”1S,L,M^M^;  ^ ?q+1)  (2^^;  x^ 


X*2  (W2  2’  V d V“ d Xn+1  * 


,.2  v , 2 


With  the  functions  , satisfying  (11)  and  with  E’  - E=  -fk^-f  kp, 
H(l,2)  + E ' - E in  (10)  may  be  replaced  by 
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(Z^-Z-n)  r"1  + (Zg-Z-n)  r2_1 


w 


V2(r2)  + r12 


-1 


(2C) 


where  r^0  is  the  interelectronic  distance. 

If  a is  orthogonal  to  co,  and  if  the  wave  functions  of  bound  elec- 
o 1 

trons  are  orthogonal  to  cDg,  which  is  always  satisfied  for  and 

l d jig,  the  first  four  terms  in  (2C)  give  no  contributions  to  the  integral 
in  (IQ)  and  can  be  omitted.  If  l = l g,  the  orthogonality  of  P^  and 
F,  , can  be  achieved  by  various  methods,  e.g.  by  redefining  the  potential 

Fk  l dr’  °r  by  solvinS 
eq.  (11)  with  the  right-hand  side  replaced  by  \ P , where  \ is  to  be 

Xj 

adjusted  so  that  I P F dr  = 0.  In  the  case  of  1=1  we  can  use 
j a kg*  l o 

similar  methods,  or  we  can  put  Z^  = Z and  V^(r)  = V (r)  in  (7)  and  (11). 

Then  cd  and  * will  be  orthogonal  because  k d k.  . In  the  following  we 
ol  o 1 

shall  assume  that  the  orthogonality  properties  are  always  satisfied  and 
that  (20)  may  be  replaced  by  r^g1. 

The  one-electron  functions  defined  by  (IS)  represent  states  of 
definite  orbital  and  spin  angular  momentum.  Using  vector-coupling 
coefficients  we  therefore  define 


V0  in  (11)  or  by  replacing 


V by 


F - P 

tkg2 


I 


•TUAl)^;  STLT1§£:  k0; 


(21) 


S i ST  1L1 


= 2-  Lm  " 'i'^nSLMCHr;X0--*XnJ)cC^nmn^k''X1)- 

MS^oMs  ^ MLmoML  Si  2 n+1  oooo  1 


sVm^  specify  spin 
follows  that 


and  angular  momenta  of  the  total  system. 


It 
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i T T 

S f S>v  L l L 


Y(  tnSL'A  ':l  ; x x )■:  ( !m~k',x  ) = ^ 

2 "+1  V.A  W£ 


(22) 


XY[(Al))o;  sVm^:  k0i  ?!— -V,i 


Similarly  to  (22)  we  can  write  in  (19)  the  product  Y(£n  ^S'L'MgM^; 


*3— *n+1)  V2  ( ^2m2^2^2 ’ V in  termS  of  Y[(An’1S'L')  Z2;  SFLFM^m£;  k2; 

x_— x ; x 1 where  S flW  characterize  the  final  state  of  the 

3 n+1  2 - * S L 


ion  + ejected  electron,  and  then  the  product  'f[(£n  1S'L')22; 
SFLFM^M^;  k2;  x ^n+l’  *2"  *1  in  terms  of 


n-1  FF  TTTT  ~+ 

Y[[(ji  S’L’)  l2\  S L*]^;  SL  MS>1£  ; k^  ; ^ *n+1 ; *2;  ^].  (23) 


The  integral  in  (10)  can  be  expanded  in  a series  of  terms  and  only  those 

TTTT 

terms  for  which  the  quantum  numbers  S L defined  by  (21)  and  (22) 

are  identical  to  specified  by  (23)  will  give  non-vanishing 

contributions . 

We  also  make  use  of  the  expansion  for  the  inverse  interelectronic 
distance, 


r"2  = kn  Z (2\  + 1) 

Xu 


-1 


X.  4*  1 


‘Xu 


<*1>  Yxu  ^ 


(2k) 


> 

and  substitute  f from  (19)  with  all  changes  into  (3).  The  integration 

A A A 

over  kQ,  k^,  and  k2  can  be  easily  performed  using  orthogonal  properties 
of  spherical  harmonics  and  the  summation  over  angular  momenta  and  their 


.11- 


components  as  well  as  spin  orientations  can  be  carried  out  using  pro- 
perties of  vector-coupling  and  Racah  coefficients.  Together  with  the 
summation  over  S'L'  we  obtain 


Z a UnSL;  k -*  in-1S'L';  k.k.)  = 1.6.17.2... 

S'L’  0 12  k0  kl  k2 

(25) 

_ , (l  L\\2(l  l \\2  2 

\ *(2X+1)“  0/  Vo  Oo)  V*Wl>  • 

ViV 


The  radial  integrals  are  defined  by 


V'Wi’  = 


/ 


W 


?k  £ <rl> 
o o 


r X 

< 


\+l  *2*2 


ko^o(r2)Fk. 


> 


A(ri> 


drldr2* 


(26) 


The  total  ionization  cross  sections  can  be  calculated  from  (1),  (2) 

and  (25).  In  (25),  the  3-j  symbols  have  non-vanishing  values  only  if 

triangular  conditions  for  vectors  l 2, X and  22_\  are  satisfied.  In 

o 1 2 

particular  that  means  that  g £ + \.  The  value  of  usually  rapidly 
decreases  with  increasing  \ and  therefore  only  few  low  values  of  the 
angular  momenta  of  the  ejected  electron  need  to  be  taken  into  account. 
The  situation  is  illustrated  in  the  next  section  for  the  ionization  of 
hydrogen. 

The  ionization  of  a complex  atom  with  n equivalent  electrons  in  the 
outer  shell  can  be  treated  similarly,  if  we  redefine  the  potentials  V^(r) 
in  (7),  (11),  (1L)  and  (16)  by  adding  contributions  arising  from 
additional  nuclear  charge  and  screening  of  inner  electrons. 


IV.  Results  of  Preliminary  Calculations 

In  our  preliminary  calculations  we  have  computed  ionization  cross 
sections  for  H,  He,  He+  and  Mg"’"  in  a distorted-wave  approximation  in 
order  to  compare  the  results  with  other  methods. 

We  have  put  = Z,  i.e.  the  asymptotic  charge  Z^  for  the  scat- 
tered electron  equal  to  the  charge  seen  by  the  incoming  electron,  and 
also  VQ(r)  = V^(r)  and  Vq  defined  by  (14),  so  that  the  orthogonality  of 
and  in  (IQ)  is  satisfied.  The  asymptotic  charge  Z^  acting  on  the 
ejected  electron  was  taken  to  be  equal  to  Z + 1 and  V (r)  given  by  (1 6) 
with  $ = 1.  The  orthogonality  of  P and  F was  achieved  by  solving 

eq.  (11)  with  the  right-hand  side  replaced  by  \P  and  with  \ appropriately 

AJ 

adjusted . 

The  results  of  the  present  calculations  are  compared  with  other 
methods  and  with  experimental  values  on  Fig.  1-5-  (Crosses  represent 
results  of  the  present  DW  calculations.  EQ  = energy  of  the  incident  elec- 
trons, E^  = ionization  energy).  Our  distorted-wave  method  is  very  close 
to  the  Born  (b)  approximation1-.  In  fact,  the  two  methods  differ  for  ioni- 

f 

zation  of  neutral  atoms  just  by  the  form  of  functions  F,  which  are  solu- 
tions  of  eq.  (11).  In  the  Born  (b)  and  the  DW  approximations  the  upper 
limit  of  integration  according  to  (2)  is  E/2,  while  in  the  Born  (a)  method 
the  limit  is  E and  therefore  this  approximation  always  gives  higher  re- 
sults for  cross  sections.  Near  threshold,  however,  the  DW  results  are 
higher  than  the  Born  (b)  values  (Figs.  1 and  h),  probably  because  of  the 
focusing  effect  of  the  short-range  nuclear  potential. 

It  is  not  clear  why  there  is  a discrepancy  between  the  DW  approxi- 
mation and  the  Coulomb-Born  tiethod  for  Mg-  (Fig.  5).  The  CB  calculation1' 


-13- 


is  based  on  atomic  wave  functions  generated  in  a scaled  Thomas-Fermi 


potential.  Moreover,  contributions  from  excitations  to  autoionizing 
levels  were  also  taken  into  account.  In  our  DW  calculations,  Hartree-Fock 
orbitals  of  all  atomic  electrons  were  employed. 

In  almost  all  cases  where  a comparison  with  experimental  results  is 
possible,  theoretical  ionization  cross  sections  are  larger  than  the 
measured  values.^  With  respect  to  many  simplifying  assumptions  made 
during  the  calculations,  it  is  not  obvious  what  is  the  main  cause  of  the 
discrepancy.  One  possible  reason  may  be  the  incorrect  choice  of  asymp- 
totic charges  and  (see  below). 

For  the  ionization  of  hydrogen,  the  total  cross  section  Q may  be 
written  as 


Q(ls,ko)  = -JL 
k 

o 


E/2 

/ £ <lsko 

•i  11-  0 


y<2)  d ($  k2  ), 


(27) 


where  Q (lsk  -*  £ k ) is  the  partial  collision  strength  corresponding 

■4  _ O c.  c. 

0 

to  incoming  electrons  with  angular  momentum  l and  to  ejected  electrons 
with  Zg  and  momentum  kp. 

In  Table  1,  the  importance  of  individual  partial  contributions  is 

shown  for  E / E.  *4  and  k_  = 0.  The  largest  contribution  comes  from 
oi  2 

^2  * 1,  i.e.  from  the  dipole  transition  s -*  p,  and  contributions  from 
higher  l 2 rapidly  decrease. 

Table  2 presents  the  dependence  of  Q (lskQ  "*  ^2k ) = T.  0, 

2 ^ a ^ 

(lskQ  -*  ^2^2^  °n  c'ie  enerSy  t*ie  ejected  electron  % k7.  For  low 

values  of  1-,  this  quantity  decreases  with  k2  indicating  that  the  most 
important  region  for  the  calculation  of  cross  sections  corresponds  to 
low  energies  of  the  ejected  electrons. 

-14- 
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The  effect  of  the  choice  of  the  asymptotic  charge  Zp  acting  on 

the  ejected  electron  is  demonstrated  in  Table  3 for  ionization  of 

hydrogen.  The  calculation  has  been  made  for  E /E.  =4-  and  with 

01 

Z = = 0.  The  Table  shows  that  a slight  increase  of  Zp  leads  to  a 

substantial  improvement  of  the  result  in  comparison  with  experiment. 
However,  at  the  present  time  it  is  not  clear  why  a specific  choice  of 
Zp  should  be  preferable  from  the  theoretical  point  of  view. 

V.  Summary 

A distorted  wave  treatment  has  been  developed  to  describe  the 
electron  impact  ionization  on  atoms  and  ions.  Preliminary  calculations 
were  carried  out  for  the  ground  state  ionization  cross  sections  of 
H,  He+,  He,  and  Mg  as  well  as  the  2s  excited  level  of  H and  are  found 
to  be  in  a good  agreement  with  the  Born  (b)  approximation^  except  at 
the  threshold  region  where  the  DW  method  yields  higher  values.  In  the 
case  of  H,  individual  and  total  angular  momenta  contributions  to  the 
partial  collision  strength  are  presented.  In  terms  of  the  individual 
angular  momentum  contributions,  at  threshold  for  the  ejected  electron, 
the  dipole  contribution  to  the  partial  collision  strength  is  generally 
the  largest,  but  not  always  as  evidenced  by  the  combination  of  ( JLq,  l„) , 
(0,1)  and  (1,0).  As  the  ejected  electron  energy  is  increased  the  dipole 
partial  collision  strength  remains  the  dominant  term  but  steadily 
decreases  in  value. 
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Table  1 

PARTIAL  COLLISION  STRENGTH  Q.  (lsk  -♦  £„k0) 

i O <L  d 

o 

FOR  IONIZATION  OF  HYDROGEN,  f = 4 ryd,  k£  = 0. 
(Numbers  in  parentheses  are  decadic  exponents) 


l2 

mm 

0 

1 

2 

3 

4 

0 

3.75C-1) 

2.03C-1) 

5.68C-3) 

3.47C-4) 

1.21C-6) 

1 

7.7K-U 

2.13 (-1) 

3-54C- 2) 

8.43C-4) 

7.89C-8) 

2 

14-.36C-1) 

1.33 

5.78(-2) 

2.12C-3) 

9 .36 (-8) 

3 

1.62C-1) 

2.01 

5.46(-2) 

3.914(-3) 

1.3CC-5) 

4 

i4-.47C-2) 

2.09 

1.12C-1) 

4.44(-3) 

2.5CC-5) 

5 

3.71 (-3) 

1.82 

1.73C-D 

/— *\ 

1 

O 

0 

7 • 33  C - 5 ) 

6 

8.99  C“i4-) 

1.44 

2.09C-1) 

4.l6(-3) 

1. 10(-4) 

7 

7A9(-7) 

1.08 

2.14(-1) 

4.93C-3) 

1.16C-4) 

8 

7.89C-D 

1 .96(-l) 

5.33C-3) 

1.10C-4) 

9 

5.68C-1) 

1 .65  C-l ) 

6.39C-3) 

1.02C-4) 

10 

4.09(-l) 

1.32C-1) 

6.92C-3) 

9.33C-5) 
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Table  2 


PARTIAL  COLLISION  STRENGTHS  (lsk  -♦  k ) FOR 

O d 2 

, 2 

IONIZATION  OF  HYDROGEN  AS  FUNCTION  OF  f k2 
(ENERGY  OF  THE  EJECTED  ELECTRON ( f k^  = 1*.  ryd  ) 


t k2  (ryd.) 

0 

1 

2 

3 

k 

0.0 

1.30 

15.01 

1.88 

C.C55 

0.0007 

0.5 

0.69 

3.57 

1.22 

0.26 

0.C26 

1.0 

0.35 

iM 

0.69 

0.2k 

O.C55 

1.5 

0.20 

0.70 

0.40 

0.18 

0.063 

I 
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Table  3 

DEPENDENCE  OF  THE  IONIZATION  CROSS  SECTION  FOR 
HYDROGEN  ON  THE  ASYMPTOTIC  CHARGE  Zg 
(FOR  THE  EJECTED  ELECTRON).  Eq/E  = 4. 


Z2 

Q Cir  %) 

0.4 

1.3^ 

0.6 

1.36 

0.8 

1.21 

1.0 

1.01 

1.2 

o.8o 

1.1+ 

0.62 

l .6 

0.1+8 

7 

Experiment 

0.69 

19 


close-coupling 
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Figure 
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